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Abstract 

We obtain a new representation for the partition function of the six vertex model with do- 
main waU boundaries using a functional equation recently derived by the author. This new 
representation is given in terms of a sum over the permutation group where the partial homo- 
geneous limit can be taken trivially. We also show by construction that this partition function 
satisfies a linear partial differential equation. 
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1 Introduction 

In the past decades the mathematical structure underlying integrable systems has been intensively 
studied and even after many years it seems that its richness has not yet been completely exhausted. 
Integrable systems can be solved by a variety of methods, ranging from functional to algebraic 
methods, whose solution is usually expressed in terms of solutions of Bethe ansatz equations [H [21 
El HI- By way of contrast, the six vertex model with domain wall boundary conditions (DWBC) is 
one exception for this general behaviour, and in fact the computation of its partition function and 
free-energy does not rely in Bethe ansatz like solutions [3 El [7] . 

In the pioneer work [5] it was demonstrated that this model obeys a recurrence relation 
relating the partition function of the system on a square lattice of size L x L to the one on a lattice 
of size (L — 1) X (L — 1). Five years later Izergin proposed a determinant solution for this recurrence 
relation which, together with extra properties, determines uniquely the partition function of the 
system [6]. 

On the other hand it is nowadays well understood the fundamental role of the Yang-Baxter 
equation and the Yang-Baxter algebra in the construction and solution of integrable systems, though 
the Yang-Baxter algebra does not seem to play any explicit role in Izergin-Korepin solution of the 
six vertex model with DWBC. Recently it was proposed an alternative approach for computing this 
partition function where the Yang-Baxter algebra is the main ingredient [8j. The approach of 
makes explicit use of the Yang-Baxter algebra in order to derive a functional equation determining 
the partition function of the six vertex model with DWBC. Though without a rigorous proof, in 
the framework of [8j the partition function is uniquely determined by three conditions: 
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(i) Functional equation 

(ii) Polynomial structure 

(iii) Asymptotic behaviour 

and here we aim to demonstrate that Korepin's recurrence relation can be suitably introduced in 
this framework removing the need of polynomial solutions. This is of particular interest for further 
applications since there exist relevant models, such as the Hubbard model, whose ii-matrix contains 
non-polynomial elements [9]. Moreover, the introduction of the recurrence relation in the functional 
equation derived in |B] yields naturally an explicit representation for the partition function of the 
six vertex model with DWBC whose homogeneous limit for the vertical degrees of freedom can 
be obtained trivially. In a second analysis we also complement the results of [8j by showing that 
the functional equation previously obtained for the partition function of the six vertex model with 
DWBC can be converted into a linear partial differential equation. 

This paper is organized as follows. In the section 2 we recall some basic definitions and 
results of [8j as a starting point. In the section 3 we demonstrate how Korepin's recurrence relation 
can fit in our functional equation yielding a novel representation for the partition function of the 
six vertex model with DWBC. In the section 4 we recast the functional equation for the partition 
function as an operator equation followed by its differential representation. Final comments and 
concluding remarks are discussed in the section 5. 



2 Functional relations and domain wall boundaries 

In this section we recall some previous results and definitions associated with the derivation of a 
functional equation determining the partition function of the six vertex model with domain wall 
boundaries. Let us consider the matrix C{X) given by 
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a{X)j 



(1) 



containing the statistical weights of the six vertex model. The non-null entries are given by a(A) = 
sinh (A + 7), 6(A) = sinh (A) and c(A) = sinh (7), where A and 7 are complex variables. The matrix 
C{X) satisfies the Yang-Baxter equation, namely 

>Cl2(A - /i)£l3(A - v)C-2z{p. = >C23(^ - i^)^13(A - i^)Ci2iX - fJ-) (2) 
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where Cij € End {Vi(^Vj) and = C^, and in this way the monodromy matrix T defined by 

T{\ {/Ufc}) = £^i(A - /xi)/:^2(A - ^i2)■■■ - AtL) (3) 

satisfies the foUowing quadratic relation usually referred as Yang-Baxter algebra 

R{\ - v) r(A, {ilk}) ® T{u, {/xfe}) = T{v, {ilk}) ® r(A, {ilk}) R{\ - v) . (4) 

Here i?(A) = PC{X) and P denotes the standard permutation matrix. 

The monodromy matrix T(A, {/Ufc}) can be conveniently written in terms of operators ^(A, {l^k})-, 
5(A,{/Ufc}), C{\,{iik}) and D{\{iik}), i.e. 

r(A.{.a)^h''*'"" (5) 

\c{\{ixk}) D{X,{iik})J 

and as demonstrated in [5j the partition function Z of the six vertex model with domain wall 
boundaries in a lattice of size L x L can be expressed as 

L 

Z = {0\llB{X„{iik})\0) (6) 

i=i 

where the states |0) and |0) consist of the ferromagnetic states 



|0) = (^y and |0) = (^^J. (7) 

Among the three conditions given in [8] determining the partition function Z, the main one 
consist of the functional equation 

L+l 

Mj Z(Ai, . . . , Aj_i, Aj+1, . . . , \l+i) 

i=l 

+ ^ ■2^(Ao, Ai, . . . , Aj_i, Aj+i, . . . , Aj_i, Aj+i, . . . , Al+i) = (8) 

^<i<j<L+l 

where we write Z{Xi,...,Xl) for the partition function ([6]) omitting the dependence with the 
variables {nj}- In their turn the coefficients Mj and Nji are given by 

c(A» - Aq) -[^ a(Ai - Afc) a(Afc - Ap) 

= KXTx^) n - '^'^'^"^ - n _ A,) 6(A, - Ao) 
i=i fc— 1 

c(Ao - Aj) -A- X TTT a(Ao - Afc) a(AA; - Aj) 

+ KX^) n - - n 5(Ao - A.) 6(A. - A.) 
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_ c(Ao - Xj) cjXj - Aq) ajXj - Xj) ^ a{Xj - \m) ajXm - Aj) 

~ b{Xo - A,) 6(A. - Ao) b{X, - A.) 1 i "^^"^ ii 5(a^. _ A^) 6(A„ - A.) 

c(Ao - Ai) c(Aj - Aq) a{Xi - Xj) _ _ '^^j a(Aj - Xm) ajXm - Xj) 

+ b{Xo - A,) 6(A, - Ao) 6(A. - A,) f^i "^^^ 1 i 6(A. - A^) 6(A^ - A,) " 

(10) 

Here we shall also employ the variables Xi = e'^^'^'-~^^'^ in order to characterize the polynomial 
structure of Z. More precisely, the partition function Q exhibits the following polynomial structure 

Z(A„...,Ai) = ^%IlllI^ (11) 



n 



L-l 
2 



i=l 

where Z(xi, . . . ,xl) is a polynomial of degree L — 1 in each variable Xi separately. 

Besides the two conditions discussed above, the full determination of the partition function 
also makes use of the asymptotic behaviour Z{xi, . . . ,xl) ~ [L]q'^^- i^i ■ ■ ■ xl) as Xi — ?> 

CX3, where [L]q2\ denotes the g- factorial function defined as 

[L],2! = l(l+q^){l + q^ + q^)...(l + q^ + ... + <^2(i-l)) . (12) 

The partition function ([6]) also exhibits some extra properties besides the three conditions 
discussed above. For instance, Z is a symmetric function under the exchange of variables Xi Xj, 
i.e. 

Z{Xi, . . . , Xi, . . . , Xj, . . . , Xl) = Z{Xi, . . . , Xj, . . . , Xi, . . . , Xl), (13) 
which follow directly from the commutation relation 

B{X, {fik})B{u, {/ifc}) = B{u, {/ifc})S(A, {f,k}) (14) 

encoded in the relation (jH). 

We close this section remarking that the partition function ([6]) is also symmetric under the 
exchange of variables fii o fij as discussed in [5j, however this property will not be required for 
our forthcoming analysis. 



3 A new representation for Z{Xi, . . . , A^,) 

In this section we aim to demonstrate how we can combine the functional equation ^ and Korepin's 
recurrence relation [S] in order to produce a new representation for the partition function Z. As 
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a matter of fact it turns out that the functional equation and the recurrence relation together are 
enough to completely determine the partition function. 

In order to construct a representation for the partition function Z we firstly recall the recur- 
rence relation derived in [5] adjusted to the notation we are considering here. 

Since this recurrence relation has been already discussed in many works [U [5l [11] we shall 
omit its derivation and here we only need the relation 

L 

Z(Ai, . . . , AlI/^i, . . . , AiL)lAi=^i = n ^^■^3 - f^iMf^i - f^j)Z{\2, • • • , Al|^2, • • • , ^l) 

i=2 

(15) 

with Z(Ai|/ii) = c(Ai — /j^i). In the relation ([15]) we have included explicitly the dependence of Z 
with the variables {/^j} and also have assumed that Z{Xi, . . . , A„,|/Ui, . . . , denotes the partition 
function Z on a lattice of size nxn. In this way for instance, the Eq. p5|) is a first order recurrence 
relation connecting the partition function on a (L — 1) x (L — 1) lattice to the one on a L x L lattice 
for a particular value of the variable Ai. 

For ilustrative purposes we shall first consider the functional equation ([5]) with L = 2 which 
is then given by 

MiZ(A2, A3|^i,/U2) + M2Z(Ai, As 1/^1,^2) + M3Z(Ai, A2I/U1, /i2) 

+iV21-Z'(Ao, A3|/Ui,/i2) + iV31^(Ao, A2|/il,/U2) + iV32^(Ao, Ai|/ii,//2) = (16) 

where the coefficients Mj and Nji follows from @ and (|1U|) with the appropriate value of L. 
By setting Aq = /^i and A3 = ^1 — 7 we obtain 

while M3, N21, N31 and A'^32 remain finite. More precisely we have 

I 2 / N MX2-H1) a(A2 - Ai) 

N3i\xo=,,M=,.-, = ^"(/^^-^iH^^-^^) a(A,-;.,)fe(A2-Ai) 
AT \ 2 / N /X 6(Ai -/ii) a(Ai - A2) 

^32|ao=,„A3=..-. = ^"(^^-^^M^^-^^) »(A,-;.,)fe(Al-A2) 

and we omit the explicit form of N2i\xq=^-^ ^3=^(1 -7 since it will not be required. 

For this particular choice of variables Aq and A3 the Eq. (|16|) is then reduced to 
, I X a(Ai - fi2) b(A2 - a(A2 - Ai) 

^(Al, A2 /il,M2) = -7 r-7T TTTT Z (/Xl , A2 Ml > ^2 ) 

a(/"i - M2) a(A2 - Ml) o(A2 - Ai) 



a(A2-M2)b(Ai-Mi)a(Ai-A2) , , . 
aim - M2) a(Ai - Ml) o(Ai - A2) 



I ^ ^21 



M3 



(19) 



-^0=Mii-*^3=A'i-7 



and if we set A2 = fJ-i in ()19p and consider that the partition function Z is symmetric under the 
exchange of variables Aj -H- Xj , we are left with the following identity 

i N21 



^0 = ''l,>'3=Mi-7 
A2=/J1 



• - 7|(Ui,/i2) = . (20) 



Since the quantity inside the brackets in ([20]) is finite, we can conclude that Z{fii, fii—j\fj.i, fi2) = 0. 
Thus the relation ()19p simplifies to 



\ I N a(Ai - M2) fe(A2 - /^i) a(A2 - Ai) 

^(Al, A2^1,/i2) = -7 r-7T TTTT -^(/^l , A2 W , /i2) 

a{fJ-i - 1^2) a(A2 - fJ-i) o(A2 - Ai) 

a(A2 - fi2) fc(Ai - //i) a(Ai - A2) , . . . 

a(/ii - /U2) a(Ai - ^1) 6(Ai - A2) 

Notice however that the property Z(/xi,/ii — 7|/Ui,/i2) = could also have been inferred from the 
recurrence relation (jlSp . 

Now we can simply insert the recurrence relation (|15p into the relation ()19p . and by doing so 
we automatically obtain 

Z(Ai,A2|/Xl,Ai2) =i^l2 + i^21 (22) 



where 



F,, = c2a(A. - /X2)6(A, - /.i) ^|^_^^j . (23) 

This procedure can be straightforwardly extended to the case L = 3. In that case the 
functional equation ([8]) reads 

MiZ(A2, A3, A4|;Ui, /i2, /X3) + M2Z(Ai, A3, Xi\ni,H2, /Us) + M3Z(Ai, A2, A4I/X1, /X2, /i3) 
+M4Z(Ai, A2, A3|/ii,/i2,^3) + N2lZ{Xo, A3, A4|^i,/i2,/^3) + A''31 ^(Aq , A2 , A4|/ii,/i2,/U3) 
+iV4lZ(Ao, A2, X3\fIi,fJ,2, /U3) + iV32-^(Ao, Ai, A4|;Ui, /i2, /^s) 

+iV422'(Ao, Ai, A3|//i,/Li2,At3) + A'43^(Ao, Ai, A2|//1,//2,A*3) = (24) 

and now we set Aq = /ii and A4 = ;Ui — 7. By doing so we find that 

The remaining coefficients do not vanish but from the recurrence relation ()15p it is easy to see that 

2'(^i,Ati -7,A|/xi,//2,/^3) = , (26) 
and we only need to consider the terms M4, N41, N42 and A'43 in the Eq. ()24p . When Aq = /Ui and 
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A4 = /ii — 7 those terms simplify to 

mi = A^4i|ao=/.i,A4=/.i-7 
rh2 = iV42|Ao=/.i,A4=m-7 



43lAo=/^i,A4=/ii-7 



and we are left with the equation 



2 / 

a(/ii - 


/i2)a(/^2 - 


Mi)a(/ii - 


M3)a(M3 - 


Ml) 




; - /il)a(/i3 


- /ii)a(Ai 


- M2)a(Ai 


- M3) 


biXo - 


ui ) 6f A'^ — 


Ui ) of Ao — 


Ai ) afAq - 


- Ai) 


' a{\2 - 


fii) a{Xs - 


Ml) ^'(A2 - 


Ai) 6(A3 - 


-Ai) 




; - ^J'l)a{^i■i 


- Mi)a(A2 


- M2)a(A2 


-M3) 


6(Ai - 


III) 6(A3 - 


Ml) a(Ai - 


A2) a(A3 - 


-A2) 


■ a(Ai - 


Ml) a(A3 - 


Ml) &(Ai - 


A2) 6(A3 - 


-A2) 


-c'^a{n2 


; - Mi)a(/^3 


- Mi)a(A3 


- M2)a(A3 


-M3) 


. &(Ai - 


ft(A2 - 


Ml) a(Ai - 


■ A3) a(A2 - 


-A3) 


a(Ai - 


Ml) a(A2 - 


Ml) 6(Ai - 


A3) fo(A2 - 


-A3) 



(27) 



Z(Ai, A2,A3|mi,M2,M3) = - ca(A2 


- Mi)«(A3 


- Mi)a(Mi 


- M2)a(Mi 


- fJ'3) — Z{X2, 


A3 M2: 


M3 


— c a(Ai 


- Mi)a(A3 


- Mi)a(Mi 


- M2)a(Mi 


- M3) — ^(Ai, 
ma 


A3 M2: 


M3 


— c a(Ai 


- Mi)a(A2 


- Mi)a(Mi 


- M2)a(Mi 


- M3) — -^(Al: 
ms 


A2 M2: 


M3 



(28) 

where we have already considered the recurrence relation p5|) . Now we only need to substitute the 
expressions ([22]) and ([23|) in ([28]) in order to obtain the partition function Z(Ai, A2, A3IM1, M2) M3)- 
With this procedure we automatically obtain 

^(Al, A2, AsImi, M2i M3) = -^123 -I- -^132 -I- -P213 -I- -P23I -I- -P312 -I- -^321 (29) 

where 

17 3 m \ /\ \ /\ \t,/\ \uf\ N a(Aj — Aj) a(Afc — Aj) a(Afc — Aj) 

Fijk = c a(Ai-M2)a(Ai-M3)a(Aj-M3)o(Aj-Mi)o(Afe-Mi)o(Afc-M2)T7v^ TTTTT TTTTT TT 

o{Aj - Ai) 0[Ak — \) o{Ak — Aj) 

(30) 

The expressions (|22|) and (|29|) suggest that for general L we can write the partition function as a 
sum over the permutation group of rather simple elements. In fact we can extend this analysis for 
general L by considering Aq = Mi Xl^i = mi ~ 7 in the relation From the definition ^ 
we can immediately see that 

M,k , =0 j = l,...,L (31) 

and from the recurrence relation (|15p we obtain the property 

Z(;Ui,;Ui - 7, Ai, . . . , Al_2|mi, • • • ,Ml) = • (32) 
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In this way, for this particular choice of Aq and A^+i and considering the recurrence relation ()15p . 
the Eq. ([8]) reduces to 

L L L 

Z{Xi, . . . , Al) = -c^ JJa(AA: - m) a(/ii - Hk)—Z{Xi, . . . , Aj_i, A^+i, . . . , Al) (33) 

j=l fc=i k=2 

where again we omit the dependence with the variables {fJ-j}- In their turn the coefficients mi and 
rhj are given by 

L 



rriL = Ml+i\ ao=mi = (-l)^^^c^ - /^i)a(/^j - //i) 



= A^L+i.l .0^.. = (-i)^c^ n "(^^ - - n !a " o'l la " x \ ■ 

(34) 

Here we remark that a similar, though not equivalent, recurrence relation has appeared previously 
in [To]. The relation (p3]) can now be iterated using the results (f22]) and ([23|) . and by doing so we 
find 

Z(Ai,...,Al)= ^ Fi„„i^ (35) 
{ji,...,jL}e5L 

where iSl denotes the permutation group of order L and 

L L L L-l L 



= n n «(^- - /^.) n ^(^- - n n _ ^-i ■ (36) 



n=l j=l j=i n=l m>n 

j>n j<n 



Since the partition function ()35p is given by a sum over permutations, it is not difficult to 
see that the invariance of Z under the exchange of variables Aj o Xj is explicitly manifested in the 
representation given by (j35p and ()36p . On the other hand, the symmetry of Z under the exchange 
of variables /Uj o fij is not apparent, though the explicit evaluation of ()35p for small values of L 
indeed corroborates this property. 

Though it is well known that the partition function Z can be written as a determinant of a 
L X L matrix [6J, it is not clear if the relations (|35p and ()36p can be converted into a determinant. 
However, since (|35p consists of a sum over the permutation group Sl, it contains the same number 
of terms as the determinant representation [6]. Furthermore, the computation of the homogeneous 
limit fik ^ fJ- from the expressions psp and ()36p is trivial, in constrast to what happens with 
Izergin-Korepin determinant representation where the evaluation of the homogeneous limit is rather 
intrincated 
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4 From functional relations to partial differential equations 

In the course of the investigation of integrable systems many connections between previously un- 
related topics have emerged. For example, the Knizhnik-Zamolodchikov (KZ) equation is a fun- 
damental differential equation in conformal field theories and its rich mathematical structure is 
manifested in the variety of topics that the KZ and its quantized version appears |12^ [T3] . Inter- 
esting enough the KZ equation estabilishes a connection between two representation theories, one 
associated to Lie algebras and the other one associated to quantum groups [13l [H] . Though not in 
the same fashion as the KZ equation connects the representation theory of Lie algebras and quan- 
tum groups, the derivation of the functional equation ([5]) explores a connection between the highest 
weight representation theories of the Yang-Baxter algebra and the 5u(2) Lie algebra [8]. Moreover, 
in conformal field theory the KZ equation is a differential equation for the matrix coefficients of 
the product of intertwining operators for an affine Lie algebra g, while in the case of the six vertex 
model with DWBC we have a functional equation for a coefficient of the Bethe vectors. 

A priori it is not clear if there exist some relation between the mentioned functional equation 
and the KZ equation or its quantized version. In order to shed some light into possible connections, 
we aim in this section to complement the results of [8j by showing that the functional equation 
previously obtained for the partition function of the six vertex model with DWBC can be converted 
into a linear partial differential equation. 

Let / be a complex valued function f{z) € C[z] and z = (zi, Z2, ■ ■ ■ , Zn) € C". For a ^ [l,n] 
we define the operator D'j' as 

A": f{zi,...,Zi,...,Zn) ^ f{zi,...,Za,...,Zn) (37) 

which basically replaces the variable Zi with Za- 

Now we shall make use of the property ()13p and in terms of operators Df the functional 
equation ^ simply reads 

Nj,D^Df+^ + [MiDf +1 + iVL+i,i A°] + Ml+1 \ Z{Xi, . . . , Al) = (38) 

which consists of a second order equation in terms of the operator Df. Moreover, as we shall 
demonstrate the operator Df possesses a differential representation when it is restricted to the ring 
of polynomials. 

In order to proceed it is convenient to define the functions 

Mi=Y[x^' Mi and Nji=Y[x,^' N,i (39) 



such that the Eq. ([38|) becomes 



i=l 



+ Ml+i } Z{xi,...,xl) = 



(40) 



where the function Z defined by the Eq. (|lip consists of a polynomial of order L — 1 in each variable 
Xi separately. 

In this way we can restrict the action of the operator Df to the space of polynomials of order 
m, and this restriction is manifested in the property 



for k > m . 



(41) 



Consequentely the Taylor expansion of / is also truncated and convergent. 

Now considering that f = f{zi, . . . , Zn) and also taking into account the Eq. ()4ip we have 

df 



f — /(^l 5 • • • 1 ^j— 1 ) ■^a ) -^i+l ; • • • 1 ^n) ~l~ 

1 d^f 



+ 



(Zi- Za) H h 



dz. 
1 d"^f 



2! dzf 

for any i € [1,?^]- On the other hand, since a ^ [1,^] we can write 



(Zj Zq) 

(Zj Zq,) 



dz'^ 



d^fi^lj • • • ) Za, Zi^i, . . . , Z„) 



02 



(42) 



(43) 



and the Taylor expansion (|42p can be rewritten as 



f{zi, ... J Zi—i, Zi, Zi-\-i, ■ ■ ■ 1 ^n) 



E 



[zi - Zaf 



.k=0 



f{zi, . . . , Zi—i, Za, Z. 



kl dz^ 

The term inside the brackets performs the operation ()37p and we thus obtain 

(Za - Zi)'' d'^ 



i+1 ) • • • ) Zn) 



(44) 



k=0 



k\ dzi 



k ■ 



(45) 



At this stage we have already gathered all the ingredients required to convert the functional 
equation ([8]) into a partial differential equation. As discussed in the section 1, the partition function 
([6]) is completely determined by the conditions (i), (ii) and (iii). Among the three conditions, 
the condition (iii) is the weakest one since it just determines the leading order coefficient of the 
polynomial Z. On the other hand the conditions (i) and (ii) play a major role in the determination 
of the partition function ([6]) and their combination is what allows us to express the Eq. ([8]) as a 
linear partial differential equation. 

Taking into account that the function Z{xi, . . . ,xl) is a polynomial of degree L — 1 in each 
variable Xi, we can substitute the representation ()45p with m = L — 1 in the Eq. ()40p and we are 
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left with 



N-. 



Xi)^{xL+l - Xjf 




l<i<j<Lk,l=0 



L L-1 



Qk 



i=l k=0 



■+1 



Xi)'' + NL+l,i{xo - Xi) 



k 



dx^ 



+ Ml+i) Z{xu...,xl) = Q 



(46) 



which consists of a Unear partial differential equation of order 2(L — 1). 

5 Concluding Remarks 

In this work we have derived a new representation for the partition function of the six vertex model 
with DWBC. After Izergin's proposal of the determinant representation, this partition function has 
been rewritten in different ways over the years [151 [T6l [T71 [T8l [19] and many connections between 
this partition function and the theory of polynomials have also emerged j20[ I21j. In particular, it 
was shown in [21j that this partition function consists of the Schubert polynomial. 

Here we have obtained a representation for this partition function in terms of a sum over 
the permutation group whose possible interpretation as a determinant is not clear so far. On the 
other hand, the representation ()35l [36|) has the advantage of allowing the evaluation of the partial 
homogeneous limit /u^ ^ /i in a trivial way. 

Furthermore, this new representation is a direct consequence of the functional equation ([5]) 
and Korepin's recurrence relation ()15p removing all other requirements considered in [8] and [5] 
concerning the nature of the solutions. It is also worthwhile to stress here that, though we have 
obtained a representation for the partition function Z from the functional equation ([5]) , we have not 
solved the equation ([8]) strictly speaking since we have only considered the Eq. ([8|) at the special 
points Ao = ^1 and A^+i = — 7. However, a rigorous proof of the uniqueness of the solution for 
the system of equations formed by ([SD and p5|) seems to imply that the representation given by 
([35]) and ([36]) indeed solves the equation ([8]) for general values of Aq and A^+i. 

Concerning a second analysis of the Eq. ([8]), we have also demonstrated in the section 4 that 
the requirement of polynomial solutions allows us to rewrite the functional relation ([5]) as a linear 
partial differential equation. The partition function of the six vertex model with DWBC is known 
to correspond to a KP r function [22j and also to satisfy a Toda lattice differential equation in the 
homogeneous limit [llj. In this way we hope the linear partial differential equation ([46p to help 
shedding some light into possible connections between this partition function and the theory of 
differential equations. 
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